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We linearize the Einstein equations when the metric is Bondi-Sachs, when the background is 
Schwarzschild or Minkowski, and when there is a matter source in the form of a thin shell whose 
• density varies with time and angular position. By performing an eigenfunction decomposition, 

we reduce the problem to a system of linear ordinary differential equations which we are able to 
solve. The solutions are relevant to the characteristic formulation of numerical relativity: (a) as 
^v^j . exact solutions against which computations of gravitational radiation can be compared; and (b) in 

formulating boundary conditions on the r = 2M Schwarzschild horizon. 
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. This work is motivated by the needs of numerical relativity, specifically the characteristic, or null cone, formulation 
\Q used in the PITT code 0, Q (and see also 0, 0, S ) . It should be noted that the PITT code is not the only approach to 
characteristic numerical relativity - for example, see 0, H, U El El • The PITT code exhibits long-term stability in 
evolving single black hole spacetimes El , and has also been successfully applied to the problem of nonlinear scattering 
of gravitational radiation by a black hole [l4j|. However, it has not been possible to compute gravitational radiation 
emitted in astrophysically interesting scenarios involving matter - for example, a star in close orbit around a black 
hole [f| . Variations of the method of calculating the gravitational radiation have been investigated 0] and 0, El , 
but the problem has not been resolved. An important difficulty in trying to analyze the problem is that, while methods 
such as the quadrupole formula indicate approximately what the final result of the computation should be, there is 
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no guidance concerning the intermediate quantities. Thus, in this paper we construct analytic solutions containing 
^Jy gravitational radiation and matter, against which a numerical code can be tested. 

In addition, we discover an unexpected issue concerning the boundary condition commonly imposed on a black hole 
horizon - the details are given in Sec. IV Al 

Linearization and eigenfunction decomposition are tools commonly used to obtain an analytic understanding of 
mathematically complex systems. The Einstein equations are a particularly complicated system of nonlinear partial 
differential equations, and there is an extensive literature on linearization tec hniq ues. Work on linearization within the 
characteristic formalism includes the Newonian correspondence method |l7l Il8l Il9l . which is related to linearization 
about a Minkowski background. More recently [2(J , the Teukolsky equation |22| has been used as a starting point 
to investigate gravitational radiation on a Schwarzschild background; but we do not use that formalism here because 
the coordinates are not those of the PITT code. The key point about the approach presented here is that we use the 
same coordinates and metric variables as in the code, so that comparisons between numerical and analytic solutions 
can easily be made. We also allow for the presence of a matter source. 

We seek a physical situation that is without symmetry and that radiates gravitationally, but in which the linearized 
problem has a solution in simple analytic form. The model problem should be relevant to the physical problem of a 
star in quasi-circular orbit around a black hole. We achieve these goals by considering a Schwarzschild black hole of 
mass M (with the case M — corresponding to a Minkowski background), surrounded by a thin spherical shell of 
matter whose density varies in time u as e lvu (with the static case recovered by setting v = 0), and whose angular 
dependence is that of a spherical harmonic. We find a solution for all metric variables in simple analytic form, except 
when M ^ and v ^ in which case solutions are obtained as power series. 

The calculations reported in this paper are very long, and would have been impossible to perform without the use 
of a computer algebra system, in this case Maple. 

In Sec. [H] we give a summary of the formalism concerning the Bondi-Sachs metric in numerical relativity, together 
with the formalism of spin-weighted spherical harmonics. Next, in Sec. we specify the model problem to be 
solved and work out the resulting Einstein equations, reducing the problem to coupled linear ordinary differential 
equations. Solutions are found in Sec. IIVI and presented according to the values of the parameters v and M. The 
paper concludes, in Sec. with a discussion of appropriate boundary conditions at a black hole horizon, and an 
outline of opportunities for taking this work further. 
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II. FORMALISM 
A. Bondi-Sachs metric 

The formalism for the numerical evolution of Einstein's equations, in null cone coordinates, is well known 0, 0, 0, 
For the sake of completeness, we give a summary of those aspects of the formalism that will be used here. 
We start with coordinates based upon a family of outgoing null hypersurfaces. We let u label these hypersurfaces, 
x A (A = 2,3), label the null rays and r be a surface area coordinate. In the resulting x a = [u, r, x A ) coordinates, the 
metric takes the Bondi-Sachs form I24LI25I 



ds 2 = - (e^(l + y) - r 2 h AB U A U B ^j du 2 

-2e 2f3 dudr - 2r 2 h AB U B dudx A + r 2 h AB dx A dx B , (1) 

where h AB h B c — S A and det{h AB ) = det(q AB ), with q AB a unit sphere metric. We work in stereographic coordinates 
x A = (q,p) for which the unit sphere metric is 

q AB dx A dx B = -^(dq 2 + dp 2 ), (2) 

where 

P = l + q 2 +p 2 . (3) 

We also introduce a complex dyad q A defined by 

/ = f(M), ?A = |(M) (4) 

with i = \J — 1. For an arbitrary Bondi-Sachs metric, h AB can then be represented by its dyad component 

J = h AB q A q B /2, (5) 
with the spherically symmetric case characterized by J = 0. We introduce the spin-weighted field 

U = U A q A , (6) 

as well as the (complex differential) eth operators 5 and 3 (see 26] for full details). 

Einstein's equations R a p = 87r(T Q ^ — }^g a pT) are classified as: hypersurface equations - i?n, q A R\ A , h AB R AB - 
forming a hierarchical set for /3, U and W; evolution equation q A q B R AB for J; and constraints Ro a - An evolution 
problem is normally formulated in the region of spacctime between a timelike or null worldtube T and future null 
infinity, with (free) initial data J given on u — 0. Boundary data for (3, U, W, J that satisfies the constraints, is 
required on T. In the case that the spacetime contains a single Schwarzschild black hole, T is usually taken to be the 
(past) event horizon r = 2M. 

In a general situation the various metric components do not admit simple physical interpretations, but nevertheless 
a review of some particular cases does illustrate the roles of the various coefficients. Minkowski spacetime is recovered 
by setting all the coefficients to zero, i.e. (3 — W = U = J = 0. Schwarzschild spacetime has J = U = 0, and is usually 
described with (3 = 0, W = — 2M, but can also be represented by (3 — (3 C (constant) and W — (e 2 ^ c — l)r — 2M . The 
quantities J and U are zero when there is spherical symmetry, and so can be regarded as a measure of the deviation 
from spherical symmetry: the quantities are interlinked, and in the linearized regime contain all the dynamic content 
of the gravitational field. 



B. Spin- weighted spherical harmonics 

We will be using spin- weighted spherical harmonics [53, using the formalism described in . It will prove 
convenient to use s Z^ m rather than the usual s Yim as basis functions, where 

s Z lm = -J= [ s Y lm + (-l) m s Y^_ m ] for m > 
s Z hn = -j= [{-l) m s Y lm - sYe-m] for m < 

sZeo = sYm, (7) 
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and note that [l4| uses the notation s Rim rather than the s Zi m used here; we use a different notation to avoid any 
confusion with the Ricci tensor. In the case s = 0, the s will be omitted, i.e. Zg m — oZi m . Some properties of the 
s Zg m that will be used later are 

• They are orthonormal 

• They are real, which is the crucial point - quantities like dU + (5U do not introduce m and — m mode mixing. 

• Defining 

L 2 = -£(£+ I) (8) 

then 

mZ lm = mZ lm = L 2 Z lm , m 2 Z em = (L 2 + 2)dZ em , mz tm = m 2 Z tm = L 2 {L 2 + 2)Z im . (9) 

• The effect of the 9 operator acting on Zg m is 



<5Ztm — \/ — L 2 l^£ m , S 2 ^m — \l ~ P ^)L 2 (^ + 2) 2 Z tm , (10) 

III. SIMPLIFICATION OF THE EINSTEIN EQUATIONS FOR A BONDI-SACHS METRIC 

A. The linearized Einstein equations 

Physically, we consider a matter distribution in the form of a thin, low-density, shell around either a Schwarzschild 
black hole or empty space. The shell is spherical at r =constant, but the density p can vary around the shell so the 
problem is not spherically symmetric. We regard the density and metric quantities as being small, i.e. 

p, J, f3,U,w = 0(e) where W = -2M + w, (11) 

with M the mass of the Schwarzschild black hole (and, of course, the Minkowski case is recovered simply by setting 
M = 0). All terms in the Einstein equations of order C(e 2 ) will be set to zero. That is, we perform a standard 
linearization of the Einstein equations. 

Using the above ansatz, we found the following forms for the hypersurface equations and the evolution equation 

R u : -p, r - 8ttT u (12) 
r 

q A R 1A - —Ud(3-2rBf] r +rdJ r + r 3 U rr + 4:r 2 U r )=8irq A T 1A (13) 
2r 

h AB R AB ■ (4 - 258)/3 + i(8 2 J + <3 2 J) + -^(r 4 9C/ + r 4 8C/) r - 2w r = %Tt{h AB T AB - r 2 T) (14) 

2 2r z 

(2M\ 
1 J r 2 J rr + 2r(rJ). ur = 8irq A q B T AB . (15) 

In addition, we evaluated the constraint equations; these will be needed only off the matter shell, and the formulas 
below are for the vacuum case 

Roo : ^(^{r-2M)w^ rr + mw + 2(r-2M)mp-Mr{W + W)-Ar{r-2M) =0 (16) 

i?oi : (^2rw. rr + 499/3 - (r 2 9C/ + r 2 9t/), r ^ = (17) 
(/'/?,, I : ^(LVd.'r, - 2i)w + 2/- 2 ( ;— 2Af)(4f + rU. r , ) +4r 2 U :- /-(<)<¥■- rt J l T ) f 2/-()J „ - 2r 1 f „, -4/'-.];.,, j n. 



(18) 
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B. Eigenfunction decomposition 



We first suppose that the various metric quantities can be written as 

J = J (r) cos(vu)B 2 Z^ m , U = U (r)cos(i/u)BZ im , = j3 {r) cos(vu)Z lm , w = w (r) coa(vu)Z im . (19) 

Of course, a general spacetime will not satisfy Eq. (|19|l . but it can be represented by summing over £ and m and 
integrating over r$ and v. However, for our purposes we regard i, m, tq and v as fixed. Later, we will express the 
time variation in the more usual way as e lvu , but for the time being we need the dependence to be explicitly real. The 
reason is that we are now using imaginary quantities as a representation on S 2 , and this is not related to the phase. 

Having chosen basis functions that are eigenfunctions of the relevant operators, the result of applying ansatz Eq. i|19|) 
is that we reduce Eqs. I|12|) to l|15fl to ordinary differential equations 

-/V cos{vu)Z lm = SttTh (20) 
r 

— (4#> - 2r/3 ,r + r 3 £/ ,rr + 4r 2 U , r + (2 + L 2 )rJ ) caa(yu)VZi m = 8nq A T 1A (21) 

2(2 - L 2 )(3 + L 2 [L 2 + 2) J + ^(r 4 L 2 U ), r - 2w , r ) cos(uu)Z em = ^{h AB T AB - r 2 T) (22) 



-2f3 + 2U r + r 2 U , r - 2(r - M)J , r - r 2 ( 1 - — ) J , rr + 2r*(J + rJ , r ) ) cos(vu)d 2 Z tm = 8nq A q B T 



r J 



where VP stands for 



AB- 

(23) 



sm(uu) 

* = -v— 7-4- 24 

cos(^u) 

The functions /3a(r), Uo(r), wo(r) and Jo(r) can be taken as real in the above equations. Thus, we can now replace 
cos(^u) by the real part of e wu , and ^cos(^u) by the real part of —ive %vu . Any imaginary component introduced 
into any of /3o( r )j Uo{r), u>o(r) or Jo(r) can now unambiguously be interpreted as a phase variation. 

C. The news 

Linearizing Eq. (B5) in 0], we find that the gravitational news is 

N= lim [ -~ r 2 J ru + ~B 2 u; + 3 2 £ ) (25) 

r^oo \ 2 2 J 

where w satisfies (see Eqs. (11), (31) and (Bl) of [l|) 

2{uj 2 + m\oguj) = 2 + i(8 2 J + 3 2 J). (26) 

Now applying the same type of ansatz to oj as to the other metric variables, i.e., 

u = 1 + to Q Z lm e wu (27) 
where luq = 0(e), we find that Eq. (126[l simplifies to 

(4c^o + 2L 2 uj - L 2 (L 2 + 2)J 0oo )Z em e lm = 0, (28) 

where Jooo means 

Jooo = lim J (r). (29) 

r — >oo 

Solving Eq. (|2"5|) for uio, we find 

L 2 

= -^-hoo (30) 

so that Eq. 1251) becomes 

N = (~ lim (r 2 J , r ) + ^ J 0oo + ft, ) vM^ " P + 2) 2 Z tm e wu . (31) 
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D. Solution procedure 

Eqs. H20|) to H23[) may be solved in the following way. First, Eq. I|20|) shows that (3q takes a constant value, say 
/3q+, in r > ro, and a different constant value, say /3o-, in r < ro (although we will set /3q- = 0). Then Eqs. H21|) 
and ill'-il) constitute a coupled system of ordinary differential equations in Uo(r) and Jo(r), and once they have been 
solved Eq. (|22[1 can be solved for wo(r). The computer algebra does solve the systems for general £, but some of the 
expressions get very complicated and in order to simplify the presentation we now specialize to the case 

£ = 2 (32) 

so that L2 = —6. In r > ro, we remove all common factors from Eqs. ()21fl and l|23(l and they simplify to (with the 
form in r < ro obtained by replacing /3o+ by zero, as well as Jq+ by Jo_ etc.) 

4/3 0+ + r 3 U 0+ , rr + ^r 2 U 0+: r - 4r J 0+ = (33) 

~2(3 0+ + 2U 0+ r + r 2 U Q+ , r - 2(r - M)J 0+<r - r 2 ^1 - J 0+ . rr + 2riv(J 0+ + rJ 0+ . r ) = (34) 

16/3 0+ + 24J 0+ - ^(r 4 f/o+Xr - 2 Wo +, r = 0. (35) 



The constraints simplify to 



^00: p(( l J -2.yi|ir„.,r-e«' l ,. -12(r - 21/ W.,+ + 12.U/-f* () - 



-4r(r - 2M)iv/3 0+ + l2r 6 ivU + + 2rivw 0+ j = (36) 
R01 ■ ^ ( rw o+,rr - 12/3 + + Gr 2 U 0+ ,r + 12rt/ +) = (37) 
q A R 0A ■ (rw 0+l r - w + + 4r 3 £/ +,r + r A U 0+ , r r + 2r 2 U + - 2Mr 3 U 0+ ,rr ~ 8Mr 2 U +. r 



2r 2 



-4r 2 iuJ 0+ - r 4 ii/U Q+ , r - 2r 2 iv(3 Q+ = 0. (38) 



IV. SOLUTIONS 



Eqs. (|33|l and l|34|l can be manipulated to give 

d 7 d 2 J 

- 2J 2 (2x + Mix 2 + w) + 2— (2x 2 + ivx - 7x 3 Al) + x 3 (1 - 2xM) — ^ = (39) 

dx dx A 

where J2{x) = d 2 Jo+ / dx 2 and x = 1/r. Clearly, two solutions for Jo+(r) are of the form C\ + C'2/r, and the remaining 
two solutions depend on the values of v and M. Once a general solution of Eq. (|39|l has been found, it is straightforward 
to calculate Uo+(r) and Wo+(r). We will present the solutions in terms of the four cases v =, ^ 0; M =, 7^ 0. 

A. Case M = v = — Static shell on a Minkowski background 



J 0+ = Ci+ + C 2+ r 



2 . C3+ . Can 



r 



tt 2 /3 + . 0(0 2C 3+ 3C 4+ 
t/o+ = h26 2 +rH ^ 



w 0+ = -lOft+r + 12C 1+ r - 6C 2+ r 3 - ^± + C 5+ (40) 



where the Cj+, i = 1 • ■ - 5, are constants of integration. An equivalent solution, with constants of integration Ci—, 
i = 1 • • • 5 and with /3q+ replaced by zero, is applicable in the region r < rg. When we substitute the above solution 
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into the constraints, we find that i?oi = is satisfied identically, and that, after removing common factors, the other 
constraints reduce to 

Roo : 4r(2/? 0+ - 3Ci+) - C b+ = 0, q A R QA : C 5+ = 0. (41) 

We require that the solution be regular at null infinity, which means that J must be bounded as r — > oo, so that 

C 2+ = 0. (42) 

The constraints lead to C\ + — |/3o+j an d since C\ + is just Jooo, this leads directly to the news TV = 0, as expected. 
In the case that the spacetime includes the origin, the metric must be regular at r = so that 

C 3 - = C 4 _ = 0. (43) 

We set the jump in /? ar r — r to be b c , and we allow arbitrary jumps in wo and Jo, r - The jump constants, together 
with the remaining unknown Cj±, are fixed by requiring continuity in the metric variables as well as satisfaction of 
the constraints. We find that the jump conditions at r — r are 

/3 + = b c , Uo+ = Uo-, Uo+.r = Uo-, r ^ Yi wo+ = wo- - 2r 6 c , Jo+ = Jo-, Jo+,r = Jo-, r , (44) 

r o 

with the jump in Uo due to the term <3/3 r in the q A R\A equation, and not sourced by the matter. We have now 
constructed a metric that is Ricci flat everywhere except on the shell r = r , and we find the matter source by 
constructing the Einstein tensor on the shell and using T a/3 = G af3 /(8tt). Every component of T a/3 is zero except 

T m = S(r-r )Z em b c /(2irr ). (45) 

If we now set b c = 2itr a pQ 1 we can interpret T af} as being the stress energy tensor of a dust shell with energy density 
p = 5(r-r )Z em p . 

B. Case v — 0, M =fc — Static shell on a Schwarzschild background 



Jo+ = 2(2M 2 - r A )C 1+ + 2 [ 2M 2 - 2vM - 2r 2 + (2M 2 - r 2 ) In ( 1 - —X) C 2+ + C 3+ + ^± + ^± 



r J J r 2 

U 0+ = 2 -^±-4(r-2M)C 1+ -4^±(2r 2 -2M 2 -2rM + r(r-2M)ln(l-^ ■ ^ + M) 



w 0+ = -4(3 0+ r + l2C 1+ r(r - 2M) 2 + l2C 2 +r 2(r 2 + M 2 - 3rM) + (r - 2M) 2 In 



+12C3+ r+ 1 -^± + C 5+ (46) 
r 

When we substitute the above solution into the constraints, we find that i?oi = is satisfied identically, and that, 
after removing common factors, the other constraints reduce to 

i?oo : 2r(/3 0+ - 6C 3+ ) + 32M 3 C 2+ + 8M/3 + - C 5+ = 0, q A R 0A ■ 32M 3 C 2 + + 8M/3 + - C 5+ = 0. (47) 

The asymptotic form of J + is 

J 0+ (r) x -2(C 1+ + 2C 2+ )r 2 + 4M 2 (C\+ + 2C 2+ ) + C 3+ + (48) 
so that regularity at infinity will require 

C 1+ + 2C 2+ = 0, (49) 

and thus J 0o o = C3+ + A)+/2- Using the constraints, it follows that C 3+ = /?o+/6, so that J 0oo = 2/3 +/3, and this 
immediately leads to the news N being zero as expected. 
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The jump conditions and stress-energy tensor at r — ro were determined as in Sec. IIV Al The jump conditions are 

2b c 

A)+ = b c , Uq+ = Uq-, Uo+,r = Uq- r H Ti 

w Q+ = wo- - 2b c (r Q - 2M), J 0+ = J _, J +,r = Jo-,r ; ttttv ( 50 ) 



r (r -2M)' 



The only non-zero components of T a $ are 



nn ( 2M\ 1 , AR M ( 2M\ 1 , R M / 2M\ 1 , . 

T» = p(l-_J ( ^T- = ^(l-— J , ^BT^ = 5^(l--j , (51) 

where we have set 

/ 2M\~ l 

b c = 2wrp 1 and p = <5(r - r )Z 2 oPo, (52) 

V ID J 

so that T 00 = p(v ) 2 and thus p may be interpreted as an energy density. The above stress-energy tensor is not that 
of a fluid, but represents a rigid shell in which the stresses are such that they preserve the equilibrium of the structure 
and prevent it from falling towards the black hole. (The above form, in which the density varies with the angle, may 
be unfamiliar, but the spherically symmetric situation is more commonplace. In this case, p in Eq. (|52|l is replaced 
by p = S(r — r )p , T 0Q and HabT ab take the forms given in Eq. l(5TJl . and qAqBT AB = 0; such a stress-energy tensor 
describes common objects like a balloon, or a bathysphere immersed in the sea.) Returning to the form in Eqs. Ij51|l 
and l|52|l . we note that the stress-energy tensor has been constructed from a metric, and therefore the Bianchi identities 
imply T°p — 0; we have checked explicitly that this identity is indeed satisfied. Since Z 2 q is negative in places, the 
stress-energy tensor Eq. <|51|) violates standard energy conditions j2^. However, that difficulty is easily rectified by 
supposing that the matter source is a sum of Zqo and Z 2 o modes, with the coefficient of Z^o larger by a factor of \/5 
than that of Z 2 o- Then, provided r > 2M, all the standard energy conditions are satisfied. 

There are 10 arbitrary constants available, and these must be chosen so as to satisfy the following 9 conditions 

• Regularity at the horizon r = 2M imposes 1 condition, C 2 - = 0; 

• Applying the constraints in the exterior imposes the 2 conditions in Eqs. (|47J) - because the jump conditions 
preserve the constraints, no additional information is obtained by imposing the constraints in the interior; 

• Regularity at infinity imposes 1 condition, Eq. (|49|) : 

• There are 5 jump conditions at r — r . 

Thus, we construct a solution that contains 1 arbitrary constant. The values of the various metric quantities at 
r = 2M are 

Uo- = fg,, U -, r = + ^, -0- = 6C0, = § + * (53) 

where Cq is the arbitrary constant, and where 

$ = ^(-3r 3 + 9r%M ~ \2r M 2 + 6M 3 ) In ^1 - — j - 14M 3 + 12r M 2 - 6r 2 M^ . (54) 

Thus, we may impose that both of Uq- and Wq-, or one of Uo- tT or Jo-, vanish at the horizon, but we cannot require 

U = w = U r = J = at r = 2M (55) 

C. Case v 7^ 0, M — — Dynamic shell on a Minkowski background 



C 1+ C 2+ e C 3+ , 2 2 . iC i+ 

Jo+ = — : tttt H 7^ [r v + 2irv — 1H 

4r 12r 3 4r 3 v ' v 

U °+ = 2^ + y^A ( 3 + 4 ^ r ) + 4r 4 ( 3 - 2lVr ) + C 4+ + — 

wo+ = ^ (1 + 2ivr) + ^-2± + 6rC 4+ ( - - r ) + Cfe+ - 10/3 +r (56) 



1/ 
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When we substitute the above solution into the constraints, we find that i?oi = is satisfied identically, and that, 
after removing common factors, the other constraints reduce to 

Ci+ - 3i/ 3 C 2 + - mrC i+ + {irv 2 - 3t/)C 5+ = 0, q A R 0A ■ C 5 + + v 2 C 2 + = 0. (57) 

Regularity at null infinity leads to 

C 3+ = 0. (58) 

Imposing the condition that the spacetime should include the origin means that a series expansion about r = of 
Eqs. H56[) must be regular. This leads to 

Ci_ = I / 2 C 3 _ ! C 2 _ = -3C 3 _. (59) 

The jump conditions at r = r<j are 

2b c Aivb c 

P0+ = O c , Uq + = Uq-, Uq + r = (To- r H s ^~ 2~ ' 

w + = w --26 c r , J 0+ = Jo-, Jo+,r = Jo-,r + r ° g C - (60) 
The only non-zero components of T"' 3 are 

T oo , ^oa ^ab ^ 2 3 2 / , , M(r-r )Z 20 e w " 

T =/, g A T u = — — , q A q B T = — — , where / = . (61) 

6 12 27rro 

In this case, as expected, the matter shell is dynamic with movement of matter within the shell (because T 0A ^ 0); 
but we do not investigate further the physical properties of the matter shell Eq. I|61|l . 

D. Case v 0, M ^ — Dynamic shell on a Schwarzschild background 

We were not able to obtain a solution in simple analytic form. However, we did obtain series solutions about the 
singular points of Eq. I|39[) . i.e. about r — 2M and r = oo. 

1. Solution about r = 2M 

Writing r = y + 2M, we found that the series solution about y = of Eq. (|39|1 contains one singular solution and 
one regular solution. We discarded the singular solution, and found 

y 2 ( UiMv-3 \ g 2 _ 

J °- = ^ SAP I 1 ~ ZM{AiM V -Z) V + ■ ■ 'J + VT2M + Cl ~ (62) 

Db- = + C 2 _ (y + 2M)3 - C 3 - (1 - - + • • • j (63) 

w - = d_6y(2 + iyv + ±iMv) - C 2 _ 3y + C 5 _ + C 3 _6(l - 2ti/M)» fl + +■■■) (64) 

y + 2714 \ 12M^ / 

The constraints were evaluated and found to be series with somewhat complicated coefficients, so at first sight it 
appeared to be impossible to satisfy them. However, the various coefficients are not linearly independent, and impose 
only two conditions on the constants Ci_, C 2 _, C 3 _ and C5-. Taking Ci_ and C 3 _ as free, we found that the 
following solution for Cs_ and C 2 _ satisfies all the constraints at least up to the order tested (i.e., 0(y 7 )) 

3(MiM 2 v 2 -8Mv + 5i + 16M 3 v 3 ) 3(8Af 3 j/ 3 - 30iM 2 v 2 - 31Mv + 6i) 

C5 --~ 63 - 2^(2 + Miu) Cl - u{2 + Miu) (65) 

„ 3i „ -38Mv - 15* + 24iM 2 v 2 + 16M 3 v 3 

C 2 - = -Cx_ - + C 3 _ — - — 66 

v bv{2 + Miv) 
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2. Solution about r = oo 



Here we use r = 1/x and the differential equation to be solved is Eq. I|39|) . This has an irregular singularity at 
x = 0, but nevertheless we were able to construct one regular series solution. Using standard techniques we then 
constructed the second solution for J 2 , but found, as might be expected, that this solution is singular at x = 0; 
therefore we discarded it. The non-singular solution is 

J Q+ = C 1+ + C 2+ x + C 3+ - ( 1 - ^-x 2 + ■■■]. (67) 
U 0+ = 2[3 0+ x - Ci+iu + C 2+ 2x 2 (l + Mx) - C 3+ x 3 ^- ( 1 - ^x + ■ ■ • J (68) 



w 0+ = -lOPo+x- 1 + Ci + 6(2x _1 + ivx~ 2 ) + C 2+ 6Mx - C 3+ 2ivx ^1 - ^-x H ^ + C 5+ . 



(69) 



The situation concerning the constraints is similar to that noted in Sec. IIVD II Here we take Ci+ and C 2 + as free, 
and find that the following solution satisfies all the constraints up to the order tested (i.e., 0(x 9 )) 

6(2. - vM) , 9(vM - 2i) 3(2 + wM) 

— P0+ 3 G 2 + T. {(1) 



3. Comments 



The number of conditions and degrees of freedom is as in the case v = 0, M ^ 0. Thus we expect that it is not 
possible to impose J = U — U. r = w = at r = 2M, unless the solution is identically zero throughout the interior; 
we have confirmed by direct calculation that this is indeed the case. 

The coefficients of the series, about both r = 2M and r = oo, have been calculated term by term, and general 
formulas for the nth coefficient are not known. While the various series do appear to have non-zero radii of convergence, 
it is not possible to make any formal statement about this issue. 

In order to work out the jump conditions and consequently the form of the stress-energy tensor on the shell, we will 
need to match the series solutions, about r — 2M and r — oo, there. Most likely this will have to be done numerically 
and is deferred to further work. 



V. CONCLUSION 



A. Discussion 



In Sec. lIVBl it was shown that, under the circumstances stated, boundary conditions Eq. (|55|l may not be imposed. 
Further, Sec. IIV D~3l indicates that this result extends to the dynamic, i.e. v ^ 0, case. Since the imposition of such 
boundary conditions is common practice 0, El EEL t ne implications of the result need further discussion. 

First, let us be clear about the conditions under which the result applies. 

1. There needs to be a matter source that is not spherically symmetric around a Schwarzschild black hole. The 
matter source must be either static, or undergoing regular oscillations (oc e lvu ). 

2. All the metric coefficients need to be either static or undergoing regular oscillations. 

Let us note here that the matter source is crucial to the argument, because a vacuum perturbation of Schwarzschild 
cannot be static or be undergoing regular oscillations: the "no-hair" theorem indicates that such a spacetime must 
relax to pure Schwarzschild (or, in general, to Kerr if the initial perturbation contains angular momentum). 

We investigate further by supposing that condition 1 above applies but not condition 2, and for simplicity we 
consider the static case. We make the gedenken experiment that a spherically symmetric shell has existed around 
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a black hole since u = — oo. At u = the matter starts to move and a new static configuration, with p a sum 
of Zoo a- n d ^20 modes, is in place from u — M until u = +00. The situation is illustrated in Fig. ^ For u < 0, 
(3 = U = U. r = w = J = at the boundary as well as everywhere inside the shell. For u > 0, changes to the 
gravitational field propagate away from the shell at the speed of light, and the metric relaxes towards a form that 
satisfies Eq. (|53[) at r = 2M. The inner boundary of the coordinate system is r = 2M on the past event horizon 
TL~ and is an ingoing null surface. Thus the effects of the shell changing cannot reach it, and the metric g a p is 
well-behaved on each null cone u =constant. However, 

lim g a p 

u — >oo 

is not continuous at r = 2M. Thus we conclude that if only condition 1 above applies, then homogeneous boundary 
conditions Eq. (|55|l may be applied for a finite time lapse it, but not for an infinite or semi-infinite time. 

In a numerical simulation using homogeneous boundary conditions Eq. I|55|l and condition 1 , the relaxation of the 
metric towards the boundary conditions Eq. Q53II . means that the metric variables must make a finite jump between, 
say, r — 2M and r = 2M + Ar. Now, Ar — > as u — > 00, so for short evolutions there is no problem with homogeneous 
boundary condions; but for longer evolutions, once Ar becomes smaller than the grid discretization, the solution will 
seem to be discontinuous and inaccuracies will develop thus making code predictions unreliable. 

The difficulty with imposing appropriate boundary data at a black hole horizon in the Bondi-Sachs formalism does 
not, to our knowledge, affect any results already in the literature. As already stated, the problem arises only if matter 
is present, and so for example it does not apply to computations involving the scattering of gravitational radiation by 
a Schwarzschild black hole |l4j. Further, the effect comes into play only as u — > 00, and so does not apply to Q where 
the computations terminate at u — O(M). The earlier work on particle evolution [f| does have long term evolutions, 
but these were not presented as physical results. 

B. Opportunities for further work 

As described in Sec. [J the results obtained here are expected to be useful in the validation and improvement of 
numerical codes using the Bondi-Sachs formalism, particularly in the computation of gravitational radiation. On this 
point, we note that the coefficient of Jooo in Eq. I|31|l scales, for large £, as £ 4 , and so it may be that noise at high 
angular frequency is being multiplied by a large number. Further investigation and numerical testing is required, the 
starting point for which would probably be the solution exterior to the shell of Sec. IIV CI (y ^ 0, M = 0). 

This paper has raised the question as to the appropriate boundary data at the horizon for a physical situation 
involving matter near a Schwarzschild black hole. However, except for the specific case of a static distribution of 
matter in the form of a thin shell, the issue has not been resolved. Clearly, a prescription is needed when there is a 
general, dynamic matter distribution. 

It is in principle possible to develop the solution obtained in Sec. IIVDI [y 7^ 0, M 7^ 0) to yield the gravitational 
radiation emitted by a star in close orbit around a black hole, because in this case the metric variables can be expressed 
as a sum of terms of the form Eq. (|19fl . While this has already been accomplished within a different formalism |30t l31| . 
an alternative derivation would nevertheless be useful. As noted in Sec. IIV CI we have not investigated in any detail 
the physical properties of a time-varying matter shell; it would be necessary to do so in order to pursue this project. 

Evaluation of the constraints (Ro a = in vacuum) has been an essential feature of this paper. Numerical codes 
based on the characteristic formalism use the constraints in setting boundary data on the inner worldtube - for example 
homogeneous boundary data (Eq. I55| l) trivially satisfy the constraints, and the constraints have been used 6], in a 
dynamic situation, to determine part of the data. In ADM type codes it is common practice to evaluate the constraints 
in the interior of the spacetime in order to monitor the accuracy of the evolution. Characteristic codes do not normally 
evaluate the constraints off the boundary (an exception being 32]), but perhaps this could lead to useful insights. 
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FIG. 1: The Penrose diagram of the spacetime described in Sec. IV Al TiT and TL + are the past and future horizons of the 
Schwarzschild black hole. The matter shell at r — ro is the straight line between I~ and I + . The line becomes dotted at u — 
indicating a change from spherical symmetry, and the thin dotted line going to TC + indicates the inner boundary of the domain 
of influence of this change. 



